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Abstract 

Several materials, such as ferromagnets, spinor Bose-Einstein condensates, and 
some topological insulators, are now believed to support knotted structures. One of 
the most successful base-models having stable knots is the Faddeev-Skyrme model 
and it is expected to be contained in some of these experimentally relevant models. 
The taxonomy of knotted topological solitons (Hopfions) of this model is known. In 
this paper we describe some aspects of the dynamics of Hopfions and show that they 
do indeed behave like particles: during scattering the Hopf charge is conserved and 
bound states are formed when the dynamics allows it. We have also investigated the 
dynamical stability of a pair of Hopfions in stacked or side-by-side configurations, 
whose theoretical stability has been recently discussed by Ward. 

1 Introduction 
1.1 Background 

Topological solitons play an important role in many areas of physics. They can be bene- 
ficial, like the Abrikosov vortices in type II superconductors or unwanted like the disloca- 
tions in nematic liquid crystals. Likewise, the creation and observation of topological soli- 
tons in Bose-Einstein condensates is nowadays routine in laboratories around the world. 
Recent experimental observations have renewed interest in the types of models that can 
support topological solitons, and it is important to understand the basic properties and 
phenomenology of topological solitons in these models. Theoretically, the possibility of 
the presence of skyrmion structures in certain materials has been known for long[T], |2] 
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(for a more recent work and references, see [3]), and this has recently been experimen- 
tally verified at least in MnSi[l] and Fei-^^CoajSi [5]. A newer theoretical discovery is 
the family of materials collectively known as topological insulators (for a review, see [6^). 
The experimental observation of topological insulator phase in Bii_a;Sba; [7j, Bi2Se3^ and 
Bi2Te3[9] confirmed these theoretical predictions and served to further add relevance to 
the basic research involving topological solitons. For the present work, perhaps the most 
interesting possibility is the discovery that frustrated magnetic materials may support 
topological insulator phases, where the wave-function in momentum-space is classified by 
the Hopf invariant ^lOj. As will be discussed below, when the system is described by the 
Hopf charge, it will support knotted structures. 

Knotted structures have a long history in physics. They were first considered by Lord 
Kelvin, who proposed in 1867 ^H] that atoms could be knotted tubes of ether. This idea 
did not yield a satisfactory atomic theory, but subsequently more realistic models have 
been proposed with potential for knotted structures. This has been done, for example, in 
the context of ferromagnets|12j. Rose- Einstein condensates [13], and optics [H]. A unifying 
feature of all these, in addition to the knotted structures, is that all these phenomena and 
their knots can be described by classical field theory. 

1.2 Hopf charge 

In a field theoretical description the knot is not made of a rope in empty space but 
through the twisting of a globally defined vector field. Here we will only consider the case 
in which the knottedness is characterized by the Hopf charge. The minimal model with 
stable knotted structures seems to be the Faddeev-Skyrme model [15] and it is believed 
that this model is contained within various models with immediate physical application, 
such as the two-component Ginzburg-Landau model modified with an additional Ward 
like term. [16] . 

The physical carrier field in the Faddeev-Skyrme model is a smooth 3D unit vector 
field = (01, 02, 03) with 0(x)-0(x) = 1. Unit vector fields with such properties have been 
proposed to exist, e.g., in super-fluid ^He in its A-phase (Continuous Unlocked Vortex) [T7] 
and in ferromagnets|12] . In order to be able to define the Hopf charge it is also necessary 
that the vector field approaches the same value at all asymptotic directions: 0(x) — )■ 0oo, 
when |x| — > oo (asymptotic triviality). Asymptotic triviality means that from the point 
of view of the Euclidean space is topologically like the sphere and this allows one 
to define the Hopf charge as an element of the homotopy class 713(5*^) = Z. A concrete 
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Hopfion vortex ring with Hopf charge Q = —1 is given by 

/'4:{2xz-y{r^ -1)) 4:{2yz + x{r^ -1)) 8{r^ ~ z^) \ 

(l + r2)2 ' (l + r2)2 ' (l + r2)2 " ) ' ^' 

where r'^ = + + z"^ , however, this configuration is not a solution for the equation of 
motion studied here. It is easy to see that = (0, 0, —1) at infinity in any direction and 
also on the z-axis, while the preimage of = (0,0,+l), which is defined as the vortex 
core, is the ring x2 + ?/2 = x, 2; = 0. 

Since the core is just a simple loop it does not explain the Hopf charge. For that 
purpose we need to consider also the preimage-curve of some other value of (j) and the 
linking of these two curves. Let us for example consider the tubular iso-surface defined by 
03 = I and on it the curve that is the preimage of = (|, 0, |), see Figure [ij The linking 
number of this curve with the core-line can now be computed using the standard rules 
of knot theory: If the top arrow can be aligned with the bottom arrow using a clockwise 
rotation, the crossing is assigned —1, for a counterclockwise rotation +1. The linking 
number is then 1/2 of the sum of signed crossings. From Figure [l] we can see that there 
are two crossings of signature —1 and therefore the Hopf charge is Q = —1. 

The Hopf charge can also be computed by integrating a charge density: For the 
antisymmetric field Fij := ■ (9^0 x dj(j) one constructs the potential Ai such that Fij = 
di Aj — djAi and then the Hopf charge is given by 

Q=^J ^'''AF.kd'x. (2) 
A proof that these two methods give the same result is given in jl8]. 

1.3 The dynamics 

For the dynamics we use a variation of the Skyrme model as proposed by Faddeev in 
1975 [I5], it is defined by the Lagrangian 

L = j [{d^4>f + 9 i^'.] d'x, F,, ■= ■ 9^0 X 9,0. (3) 

Note that in the static case we have the following scaling property: under the scaling 
r — )■ Ar the integrated kinetic term scales as A while the integrated F'^ term scales as A~^. 
Thus according to the virial theorem nontrivial configurations will attain some fixed size 
determined by the dimensional coupling constant g. In the computations reported below 
we used (7 = 1/2. 

A very important property of this model is that the energy is bounded from below [19] 
and from above [20] 

h'\Q\^ <E<h\Q\i, (4) 
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where /i, h' are some constants and Q the Hopf charge. Since the exponent of the charge 
is less than 1 it follows that configurations with higher Hopf charge tend to form bound 
states rather than split into several charge 1 states. 

An intriguing question is the form of the minimum energy states with higher Hopf 
charge. After some tentative works [2T1 [22] detailed taxonomy has been obtained by 
Battye and Sutcliffe[2Sl 1211 125] and Hietarinta and Salo[2Sl I2Z]- The results follow very 
well[271l25] the bound Q. It turns out that a twisted torus-like configuration (similar to 
([T])) is a minimum energy state only for \Q\ < 4 [27J, for higher charges linked or knotted 
configurations are found. The trefoil knot is persistently found for charge |Q| = 7 from a 
wide variety of initial states, even for long vortices. [2H] 

We will next give the computational method and then describe the results for two 
dynamical studies: scattering head on collisions, and attractive and repulsive channel 
simulations inspired by the work of Ward in |29]. 

2 Computational method 

For scattering dynamics we used the Lagrangian ([s]) in the full 3+1-dimensional Minkowski 
space-time (with c = 1) and the resulting Euler-Lagrange equations, derived with a La- 
grange multiplier to guarantee 0-0 = 1 during evolution. The second order (in time) 
equation was split into two first order equations for numerical computations. The spa- 
tial derivatives are discretized on 5'^ points. The formulae are quite similar to those in 
j3U] . Time evolution was computed using 4th order Runge-Kutta (five stage low mem- 
ory version[31j)- For scattering studies the computational lattice varied between 200^ and 
300^, but all the simulations used the same lattice constant h = 0.15. In the three-Hopfion 
collisions the physical size of the grid was 45.0 x 30.0 x 30.0 and in the double-collisions 
30.0^. The constant of Courant-Friedrichs-Lewy like condition linking time and space dis- 
cretization is Cfi = At/h = 0.02. Ward's channel simulations were computed in a lattice 
of 300^ unless otherwise noted; the grid constant was still h = 0.15 but c// = 0.2. The 
Courant instability was controlled by smoothing. For computational platform we used 
SAMRAI[32]. 

As our work progressed the significance of boundary condition became more and more 
clear and thus we have tried a variety of different methods and parameter configurations to 
minimize boundary errors. The presented two- and three-Hopfions scatterings are calcu- 
lated with Dirichlet boundary conditions 0|a = 0oo = (0,0,-1), and the Ward's channel 
simulations with Neumann boundary condition d{x^y^z}(p\d = (0,0,0). Furthermore, in 
the two- and three-Hopfion scatterings we had a thin and mild absorption layer near 
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the boundary in order to prevent radiation from reflecting and re-entering the scattering 
region. 

Visuahzation of a vector field is problematic. Since we cannot plot the vector field 
at all points of the 3D space we chose to plot iso-surfaces determined by the value of 
the third component 03, often with 03 = 0.7, which is a narrow tube around the core. 
This tube is painted[2S] with color determined by the first two components of (f). For 
visualization platform we used Visit [33] . 

3 Results 

3.1 The fundamental deformation processes 

Before analyzing specific scattering processes it is important to recall the possible elemen- 
tary deformation processes. It was already noted above that it is not enough to consider 
only the core-line but also how the unit- vector field twists around it. In fact, the proper 
knot theoretical setting is to use framed links which can be realized as directed ribbons 
(pre-images of line segments |34j). This adds local information near the curve, e.g., twist- 
ing around it. Knot deformations correspond to ribbon deformations, which allow certain 
types of crossing and breaking, in which the Hopf charge will nevertheless be conserved. 
Examples of ribbon deformations can be seen in|34j. 

3.2 Collisions and scattering 

We will now consider the scattering of low-charge "un-knot" Hopfions in two different 
situations: three-body scattering on a fixed target with projectiles from left and right, 
and two-body scattering with impact parameter. The the individual Hopfions (minimum 
energy states) were first created using the steepest descent method (corresponding to 1st 
order dissipative dynamics with a fictitious time) [26] . These Hopfions were then rotated, 
Lorentz boosted and combined to form the initial configurations for scattering. This was 
done using interpolation routines provided by the SciPy python library. 

In the first scattering process, for which snapshots are given in Figure [2| with full de- 
tails in [Movie Sl[ we have two charge Q = —2 Hopfions approaching a stationary Q = —3 
target head on, from left and right. The moving Hopfions were Lorentz boosted numer- 
ically to speeds 0.75 and —0.75, respectively, and then immersed into the computational 
box. Since the total charge is Q = — 7 we expect the result to be a trefoil knot. This is 
indeed the case, after several intermediate steps. The resulting trefoil is not stationary 
but rather oscillates and rotates, and eventually the Hopfion settles down drifting slowly. 
Since 7^/7(2 ■ 2^/4 + 3^/4) ^ 0.76 there is theoretically 24% excess energy, some of which 
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will radiate as low amplitude, close to vacuum (03 ^ —1), waves from the scattering 
region. These are not visible at the iso-surface 03 = 0.7 level. The resulting drift is 
probably due to momentum conservation in an asymmetric radiation background. 



In Figures (for details see Movies S2j and S3) we have a situation similar to the 



one in Figure |2j except that now the central target has charge Q = +3. The total charge 
is then Q = —1 and indeed the final state has one Hopfion of Q = —1. As the Hopfions 
approach each other they deform by the rule "like colors attract", thereby preparing for 
ribbon crossings. Since 1/(2 ■ 2^/4 + 3^/4) ^ 0.18 most of the energy has to radiate away. 
This can be clearly seen in Figure |4] (Movie S3) which describes in further detail the later 
stages of the process. In this figure we have indicated the distribution of energy density 
by a white cloud. Up to t = 20.1 the energy is localized near the vortex tube but when 
the deformation towards Q = —1 continues more radiation is released and can be seen to 
radiate away. 

For the two-body scattering with impact parameter we present the scattering of two 

we 



Q = —3 Hopfions with two different impact parameters. In Figure [5] (Movie S4) 
have the impact parameter 3.75 and the Hopfions have speed 0.5. Although the Hopfions 
touch briefly at t = 12 the speed and distance prevent a bound state from forming. The 
individual Hopfions just continue tumbling along a slightly bent trajectory. 

For the same speeds but with slightly smaller impact parameter value 3.6 the result is 
entirely different (Figure |6| Movie S5). For a moment it seems that the Hopfions would 
again continue along their individual trajectories but there is just enough time to form an 
elongated loop. After this the evolution is typical for the total charge |(5| = 6 case: The 
loop behaves like an over-twisted rubber band and proceeds to make one ribbon crossing 
deformation to reach the linked loop configuration that is standard for this Hopf charge. 
The linked configuration continues to rotate and oscillate. 

The above results illustrate the difficulty associated with trying to make any gener- 
alizations about the behavior of Hopfions. However, some progress has been made by 
modeling Hopfions with elastic rods[3S]. It would be interesting to see if the elastic rod 
model could be extended to approximate dynamical features as well. 

3.3 Static initial states 

We will next consider the case of two unit charge Hopfion unknots placed on static po- 
sitions not too far from each other. This case was studied by Ward [29], who observed 
(using energy minimization) that under symmetry preserving relaxation an axisymmetric 
two Hopfion state reaches a local minimum [E = 2.26) where the soliton cores around the 
symmetry axis are separated by a nonzero but finite distance. He also noticed that this 
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separated configuration has 13% higlier total energy than the twice wound multisohton 
state [E = 2.00) both of which are in the same homotopy class. The nature of the ax- 
isymmetric state is not yet clear and, as Ward stated, more study is needed to investigate 
whether the state is actually a local minimum or just a saddle-point. 

For these cases we prepared the initial states using Wards proposal, that in the stere- 
ographic coordinates W = (0i — z02)/(l + ^a), a good ansatz for a |Q| = 1 Hopfion is 
given hj W = (x + iy)/{z — if(r)), where r = (x^ + + z'^Y^'^ and / is some specific 
function of r (see Equation (4) in [29j, Ward uses W = (0i + i02)/(l + ^s)). Actually this 
can be generalized by adding some signs in the expression: 

w= "■"+'"^';, . (5) 

z - iTaia2f{r) 

Here the sign of r determines whether this is a Hopfion or anti-Hopfion, while the signs 
of the (j's can be changed by 180° rotations around one of the coordinate axes. Thus for 
a single Hopfion the signs are irrelevant but they do play a role when we construct an 
initial state with two Hopfions. 

3.3.1 Stacked Hopfions 

It was noted by Ward [29\ that a good two Hopfion ansatz is obtained just by adding the 
corresponding W functions ([s]). Thus for the case with two Hopfions located at z = ±a 
and sharing the same axis (Ward's Channel A) we can take 

^ (TjX + ia2y ^ cr[x + ia'^y 

z-a- iTcria2f{r^_a) z + a- ira[cxyXr,+a) ' 

as the starting configuration. Here rz±a = {x"^ + y'^ + (z ± a)^)^/^. 

Ward uses a dipole pair to describe the configuration, in our notation P+ = (cri, 0, 0), 
(0,^2,0) and P = (ct;,0,0), Q_ = (0,^2,0), so that := x = (0, 0, (Ti(T2), etc. 
Furthermore the angle between P+ and P_ is if ai = a[ and 180° if o"i = —cr'i- 

The signs in ^ can be changed by various 180° rotations and this way the 16 possible 
sign combinations (ai, cr2, o-[, (Xg) can be reduced to a canonical 6. They are 1) (+ + ++) 

(D's to same direction, angle 0°), 2) (+ H ) (D's to same direction, angle 180°), 3) 

(+ + H — ) (D's away from each other, angle 0°), 4) (+ H h) (-D's away from each other, 

angle 180°), 5) (H 1-+) (-D's towards each other, angle 0°), 6) ( — h ++) {D's towards 

each other, angle 180°). 



^Note that in the asymptotic hmit the situation is simpler because from the expression 

^ _ aix + ia2y _ X + icria2y 

one cannot extract r in order to say whether it describes a Hopfion or an anti-Hopfion. 



7 



The time evolution for these cases is as follows. For Cases 1 and 2 the system is in 
the attractive channel (Figure [7| Movie WIJ . The Hopfions approach each other, bounce 
a few times, during which they develop some asymmetry and eventually turn over, join 
and deform into a Q = 2 single ring Hopfion. The four remaining cases 3, 4, 5 and 
6 are repulsive: As the Hopfions recede they slowly turn around their axis and it is 
plausible that they eventually turn enough to attract each other and return to form a 
single ring state. That is, we know that due to the fact 2 x 1^/^ > 2^/^ the asymptotically 
separated configuration is a higher energy state than the single-ring configuration but it 
is unclear whether the initial configurations for these four repulsive systems have large 
enough repulsive (escape) energy to form the higher energy state. As the Hopfions quickly 
reach the border of any finite computational space, this is a hard question to answer using 
numerical methods alone. 

3.3.2 Boundary effects 

In order to better understand the effect of boundaries, we considered the (+ H ) case 

further and calculated its evolution in physical grids of different size, while keeping the 



grid constant and boundary conditions unchanged. Our results are given in Movie W2 
which combines the result from cubes of physical size 30.00^, 45.00^, 67.50^, and 101.25^ 
having 200^, 300^, 450^, and 675^ grid points, they are located at top left, top right, 
bottom left and bottom right, respectively. 

In all cases the simulation starts with a slight repulsion during which the colors of the 
isosurfaces rotate towards the same phase. Then the Hopfions approach each other and 
bounce several times and eventually form a Q = —2 un-knot Hopfion. 

The difference in the dynamics can be explained by the existence of two competing 
effects: the fixing of the boundary value to = (0,0,-1) tends to stabilize the con- 
figuration, while the radiation reflecting from the boundaries tends to un-stabilize it. 
The radiation effects are illustrated in Figure [8| where we have plotted the values of the 
usually very revealing quantity {dtcf)) on a slice of the box. The box in upper cases is 
relatively small and the closeness of the flxed boundary stabilizes the situation well and 
many bounces are required to un-stabilize it, even though there is lots of radiation bounc- 
ing around. The fastest un-stabilization occurs with box size 450^ (lower left), in that 
case radiation effects overwhelm the boundary stabilization. In the largest box (lower 
right) the radiation takes a longer time to reach the interaction region and therefore more 
bounces are again possible before deformations. 

Clearly the major difference between these simulation is the time at which the radiation 
wave of the initial relaxation travels to the boundary and is reflected back to center of the 
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grid. In smaller grids the boundary is so close that the first reflection happens before the 
Hopfions even collide the first time. After a few more reflections an interference pattern 
is formed in the background but it is reasonably symmetric. In the larger grids the first 
reflected wave passes over the Hopfion pair during the first couple of oscillations and this 
asymmetric perturbation causes the state to unwind faster. Nevertheless the final result 
is always the same 

3.3.3 Side-by-side Hopfions 

Another case considered by Ward is that of Hopfion rings side by side on the x-axis. The 
ansatz is now 

^ ^ ai{x - a) +ia2y _^ a[{x + a) + ia^y 

Here rx±a = {{x ± a)^ + y"^ + z'^Y^'^. In this case also the 16 possible sign combinations 
can be divided into 6 essentially different ones, to which others are related by rotations: 

1) (+ + ++), polarization vectors D parallel, 2) (+ H h), antiparallel, 3) (+ H ), 

parallel, 4) ( h+), parallel, 5) (+ + H — ), antiparallel, 6) (H antiparallel. 

Cases 3 and 4 (and Cases 5 and 6) are related by an overall sign change, which cannot be 
generated by a rotation. 

Considering Hopfions asymptotically as dipole pairs and studying the forces between 
them led Ward to introduce the second attractive channel which corresponds to our Case 
3 above. Ward's energy minimization gave an axially-symmetric ring for Q = 2. This is 
believed to be the true global static minimum. [29] 

Our numerical simulations show that the attractive and repulsive channels yield the 
same qualitative behavior as in the stacked case. However, now the attractive channel 
comprises of initial configurations 3, 4, 5, and 6 while the repulsive channel contains 1 
and 2. As an example of the attractive channel, we give case 5 in Figure [9| [Movie W3| . 



In order to study the repulsive case further we chose to change somewhat the rotation 
angle of the initial state. After rotating the second Hopfion by angle 9 the ansatz ([T]) can 
be written as 

ai{x — a) + ia2y a[(^{x + a) cos9 + ysm9) + ia2(ycos9 — (x + a)sm9) 
^ = ■ 77 ^ ■ — TTTT \ • (8) 

As can be seen the dependent variable rx±a is invariant under rotation. 

Recall that the configurations (+ + H — ) and (+H h) were respectively in the attrac- 
tive and repulsive channels. Now let us fix a configuration to repulsive (+ H h) with 

^ = 0, then the attractive (+ + H — ) is obtained from ([s]) with 9 = 180°. It has already 
been observed that even in the repulsive case the Hopfions turn around and could return 
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if the computational box would be big enough. In order to study a case with only slight 

repulsion we chose the configuration (+ H h) for 6 = 20° and found it to be attractive 

(Figure 11 and [Movie W4 ') while = 5° was still repulsive (Movie W5). Similarly, for 

states (+ H ) and (+ + ++): (+ + ++) is still repulsive at ^ = 3° but changes to 

attractive at = 10°. 

It seems the attractive channels yield qualitatively the same behavior regardless of 
whether the configuration is stacked or side-by-side and likewise for the repulsive case. 
The study of the dependence on the initial rotation angle, however, shows that there are 
quantitative differences, of which we have uncovered one: the repulsive channel contains 
a different range of initial relative orientations depending on the parameters a and a'. 
This further illustrates the extremely rich set of phenomena associated with Hopfions. 



4 Conclusions 

Now that we are close to seeing knotted structures in experiments it is important to 
understand their dynamics. The static knots in the Faddeev-Skyrme model are already 
well understood and in this paper we have shown that these objects have also dynamically 
many particle-like properties. 

We have here studied the scattering collisions of low charge Hopfions. In general the 
results show remarkable particle-like behavior. We have verified again their tendency to 

3 

form bound states, due to the charge-energy dependence E ~ IQI*; we have also seen the 
associated radiation of the released energy. With suitable speeds and impact parameters 
we can also have grazing scattering, without bound state formation. Since the internal 
structure of Hopfions is quite complicated the results of the scattering process depend 
strongly on the initial configuration, and not only on the initial positions and velocities 
but also on the internal specifics (such as orientation, spin, and relative phase) of the 
scattering objects [36]. 

We have also investigated the dynamical stability of static configuration of two stacked 
or side-by-side unknots discussed by Ward [29j. As expected their evolution depends 
sensitively on computational aspects due to the mere length of the calculation (due to 
the initial state being static). The errors arise both from discretizations and in particular 
from boundary effects of the finite computational box. It seems that the stacked and 
side-by-side configurations made of two charge 1 Hopfions are all unstable and eventually 
relax into a charge 2 Hopfion. 
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Figure 1: Red line: the preimage of = (0,0, +1) (vortex core); blue line: the preimage 
of = (|,0, |). The lines have been assigned compatible directions and their linking 
number can then be computed as explained in text. 
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t = 45.00 



t= 75.00 



Figure 2: The formation of a trefoil knot by the scattering of two Q = —2 Hopfions with 
velocities ±0.75 on a stationary Q = —3 target. For more details see Movie Sl[ 
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Figure 3: Collision of three Hopfions of charges Q = —2, +3, —2 and velocities 
0.75, 0, —0.75, respectively, resulting with a Hopfion of charge —1. For details see Movie| 
1521 
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t=45.0 




Figure 4: Same scattering process as in Figure [3] but with further intermediate states for 
the later stage of the deformation process. Here we have added the energy density plot 



which appears as a grayish halo. See also Movie S3 
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t= 12.00 



t= 19.80 



t = 30.00 



O 



t = 40.20 



Figure 5: Collision of two Hopfions of charge Q = —3. The initial velocities were ±0.5 



and the impact parameter 3.75. For more details see Movie S4 
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t = 0.00 



t= 12.00 



t= 19.80 



t = 37.50 




t = 49.50 



t = 63.00 



Figure 6: Same as in Figure [5] except that impact parameter is 3.6. For more details see 
fMovie S5. 
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t=330.0 t=336.0 t=387.0 



Figure 7: Stacked configuration (+ + ++) produces an attractive interaction between the 
two H = 1 Hopfions. The pair oscillates multiple times before a slight asymmetry brings 
the system to the more stable H = 2 single-ring configuration. For details see [Movie Wl| 
and the discussion in the end of section I3.3.1I 
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Figure 8: Figure shows four 30^ slices of the canonical derivative field il)2 between a fixed 
interval [—0.05, 0.05] at the same moment [t = 60.0, # = 2000) in different physical size 
grids. We can see that the boundary effects are strongly present at the area around the 
Hopfion's core in the first case, just reaching it in the second and third case and still 
outside the observed area in the fourth case. 
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t = 27. 



Figure 9: The side-by-side configuration (+ + H — ) seems to be attractive without delay 
and the pair end up in the famihar H = 2 single ring configuration. At the end of the 
simulation the pair oscillates relatively softly having already emitted most of its excess 
energy. For more detail see ,Movie W3, 
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Figure 10: The side-by-side configuration (+ H h) is repulsive from tlie start of tlie 

simulation and the system end up in the unattached state where the initial Hopfions 
advance in opposite directions. 
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t = 87. S 



Figure 11: The behavior of the side-by-side system (+ H h) changes after a rotation of 

the second Hopfion by an angle of 20°. Now the Hopfions move further away from each 
others under a color rotation until the phase difference between them becomes attractive 
and the single ring configuration if = 2 is formed. For more detail see Mo vie W4| ). 
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